Mon. Not. R. Aston. Soc. 000.mri4l(2012) Printed 1 1 December 2012 (MN LSTgX style file v2.2) 



The structure and kinematics of dense gas in NGC 2068 
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ABSTRACT 

We have carried out a survey of the NGC 2068 region in the Orion B molecular cloud using 
HARP on the JCMT, in the 13 CO and C 18 (7=3-2) and H 13 CO + (J =4-3) lines. 
We used 13 CO to map the outflows in the region, and matched them with previously defined 
SCUBA cores. We decomposed the C 18 and H 13 CO + into Gaussian clumps, finding 26 and 
8 clumps respectively. The average deconvolved radii of these clumps is 6200 ± 2000 AU and 
3600 ± 900 AU for C ls O and H 13 CO + respectively. We have also calculated virial and gas 
masses for these clumps, and hence determined how bound they are. We find that the C 18 
clumps are more bound than the H 13 CO + clumps (average gas mass to virial mass ratio of 4.9 
compared to 1.4). We measure clump internal velocity dispersions of 0.28 ± 0.02 km s _1 and 
0.27+0.04 km s _1 for C ls O and H 13 CO + respectively, although the H 13 CO + values are heavily 
weighted by a majority of the clumps being protostellar, and hence having intrinsically greater 
linewidths. We suggest that the starless clumps correspond to local turbulence minima, and 
we find that our clumps are consistent with formation by gravoturbulent fragmentation. We 
also calculate inter-clump velocity dispersions of 0.39+0.05 km s _1 and 0.28+0.08 km s _1 for 
C 18 and H 13 CO + respectively. The velocity dispersions (both internal and external) for our 
clumps match results from numerical simulations of decaying turbulence in a molecular cloud. 
However, there is still insufficient evidence to conclusively determine the type of turbulence 
and timescale of star formation, due to the small size of our sample. 

Key words: ISM: clouds - ISM: individual (NGC 2068) - stars: formation. 



1 INTRODUCTION 

An understanding of the properties of molecular gas is important in 
the investigation of star formation in our Galaxy. A complete theory 
of star formation should be able to explain certain basic properties 
and characteristics of dense cores and protostars. The kinematics 
of the dense cores and their surrounding natal environment, as well 
as the relation of the Core Mass Function (CMF) to the stellar Ini- 
tial Mass Function (IMF), are all linked to the initial condition s of 
star formation. Most stars form in clusters dLada & Ladall2003l) . al- 
though |Bjessert_eta]J yOlOj) find that the percentage of stars found 
in clusters appears to be dependent on the definitions used. Many 
studies of the less dense gas ( e.g. CO low-J transitions) i n thes e 
clusters exi st (e.g . lOka et all \ 19981) , |Pi Francesco et alj J2002h . 
iDu & Yand d2008h ). However, there are fewer surveys of the very 
dense gas tracers (e.g. HCN or HCO + 7 = 4-3 transitions), which 
reveal the density peaks in the gas most closely associated with the 
protostars. With the advent of receivers like HARP on the JCMT 
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jBuckle et alj|2009h . and HERA on IRAM dSchuster et alj|2004h . 
large-scale dense gas surveys are now feasible. 

In recent years, a general consensus has been reached that star 
formation is a dynamical process, occurring within a few freefall 
times; the method by which molecular gas clouds collapse to form 
protostellar cores is also generally accepted to be a combination 
of gravity and turbulence. This gravoturbulen t fragmentation is es- 
sentially a two-phase process jKlessenll20Tlh : the turbulence pro- 
duces density enhancements in the molecular gas cloud through 
a network of interacting sho cks; these density enha ncements de- 
crease the Jeans mass locally dClark & Bonnel3l2005l) allowing the 
densest regions (formed at the shock intersection points) to col- 
lapse, accrete mass and form stars. It should be noted that although 
these turbulent shocks are expected to be widespread in molecu- 
lar clouds, there is no observational evidence for them as yet. The 
overall picture is complicated by the presence of magnetic fields, 
which play a role on la rge scales, and cau se modifications to the 
core collapse processes jPrice & Batej|2008l) . 

Although turbulence plays such an important part in star for- 
mation, its origins, length- and timescales are still not understood. 
It is unclear whether turbulence in star formation regions decays 
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freely in a crossing time l lElmegreenll2OO0l) . or is driven dKlessenl 
1200 ll) and therefore persists f or a longer timescale . Most com- 
petiti ve accretion simulations dBonnell & Batdl2006l : ISmith et al .1 
120091) feature molecular clouds initially supported by turbulence 
that decays on a timescale comparable to the crossing time of the 
molecular cloud. This allows competitive accretion to start on small 
scales and grow progres sively larger, fe d by loss of kinetic energy 
from the ambient cloud. iKlessenl J200lh however have shown evi- 
dence for competitive accretion in an environment with driven tur- 
bulence; as long as the turbulence is driven on large enough scales, 
clusters can form where gas and stars are virialised. They have also 
produced representative core mass functions from their numerical 
simulations; however they find that with these alone they are un- 
able to distinguish between the driven and decaying modes of tur- 
bulence. Indeed, the general properties of the cores produced in 
differing turbulent environments at low resolu tion are insufficently 
different dOffner. Hansen. & Krumholzll 20091) to dismiss a particu- 
lar environment. 

We have therefore decided to use an alternative method — ex- 
amining the kinematics of the dense gas cores, rath er than their spa- 
tial and structural properties. lOffner et alj j2008j) have performed 
comparisons of the kinematic properties of protostellar and starless 
cores produced in simulations that utilise both driven and decay- 
ing turbulence modes. Comparing the intercore velocity dispersions 
(i.e. the dispersions of the core centroid velocities), they found that 
protostellar cores in decaying turbulence simulations have a higher 
intercore velocity dispersion than the ambient gas dispersion; pro- 
tostellar cores and starless cores in driven turbulence simulations 
have subvirial intercore velocity dispersions and intercore velocity 
dispersions close to that of the gas respectively. They therefore sug- 
gest that comparing the starless and protostellar intercore velocity 
dispersions to the net gas dispersion could potentially distinguish 
between the two environments observationally. 

We have chosen to observe the H 13 CO + 7=4-3 transition 
as it is optically thin, traces the densest regions (n cril ~ 10 6 cirr 3 ) 
and will be most likely to trace the dense inner regions of prestel- 
lar cores. We can calculate linewidths and virial ratios of H 13 CO + 
cores to determine how well bound they are; we can also potentially 
use the core-to-core velocity dispersions to distinguish between the 
two turbulence scenarios. We have also observed the C 18 7=3 
- 2 transition which is optically thin over most of the region and 
will also trace more of the less dense core envelopes for compari- 
son with the H 13 CO + . In addition, we observed the 13 CO 7 = 3-2 
transition to investigate the large-scale kinematics, as it should trace 
more of the bulk gas in the region than the other two molecules, as 
well as providing evidence for outflows which would pinpoint pro- 
tostellar cores. 



1.1 NGC 2068 region 

The Orion B molecular cloud complex is the closest region of 
high-mass star formation, and lies at a dista nce of around 415 pc 
dAnthonv-Twarog|ll982l ; iMenten et alj|2007t) . It includes 5 major 
star formation regions, one of which, NGC 2068, is the subject of 
this paper. NGC 2068 is a bright reflection nebula containing an 
infrared cluster; submillimetre emission associated with the nebula 
has been observed to the south, as shown in Figure [T] 

The fi rst unbiased molecular gas s urvey of NGC 2068 was car- 
ried out bv lLada. Bally. & Starkl dl99ll) . hereafter LBS90, using CS 
7 = 2-1 emission. They identified 2 clumps — LBS 10 and LBS 
17 — at 1.7' resolution. The region has since been mapped by m any 
others, both in molecular line emission (e.g. llkeda et alj J2009h . in 
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Figure 1. A DSS2-IR near-infrared image of NGC 2068, with our C ls O 
emission contours overlaid, integrated over velocity range 8.5 - 12.5 
km s . Contours are from 2 - 9 K km s , at 1 K km s intervals. 



H' 3 CO + 7 = 1 - 0) as well as i n SCU BA 850 [im dust emission 
(e.g. iNutter & Ward-Thompsorj ( I2007I) . hereafter N WT), breaking 
the tw o LBS regions into smaller clumps and cores. Mitchell et aU 
d200ll) identified 18 smaller clumps with diameters of ~ 25" using 
both 850 (im dust emission and CO isotopologue mapping. The 
clumps they identified are labelled OriBsmm (for Orion B subm- 
millimetre emissio n) and then numbered by decreasing declination. 
iMotte et al.ld200lh . hereafter M01, identified 31 cores withFWHM 
diameters of ~13" from 850 p.m dust emission; their cores are la- 
belled LBSX-MMY, where X is the corresponding LBS clump that 
the core belongs to (either 10 or 17) and Y begins at 1 and increases 
with RA. 

The NGC 2068 region has been extensively studied and well- 
documented at different densities and resolutions, and is a fairly 
small, self-contained region with obvious dust and gas structure. 
We have therefore chosen this region for investigations of the kine- 
matics of pre- and protostellar cores and their surrounding gas, in 
order to examine the initial conditions for star formation and deter- 
mine the possible mechanisms for star formation in this region. 



1.2 Outline 

We present an analysis of spectral line data from NGC 2068, ex- 
amining the gaseous structure of the region, as well as its kine- 
matics. We break the region down into Gaussian clumps using two 
molecules that trace different densities, to examine the kinemat- 
ics and mechanisms leading to the fragmentation and formation 
of prestellar cores. Section [2] gives an overview of the observa- 
tions and data reduction procedure. Section[3]describes the reduced 
data sets, highlighting the multiple-velocity-component structure of 
the gas, and several high-velocity outflows. Section|4]presents the 
clumpfinding analysis of the region, using C ls O and H 13 CO + , while 
Section|5]examines the characteristic properties of the clumps, in- 
cluding shape and size, masses and velocity dispersions. We also 
present a comparison of observational data with numerical simula- 
tions of star formation. We summarise the results and conclude in 
Section|6] 
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Table 1, Molecules observed in this particular run for NGC 2068, with their 
transitions (Column 2), frequencies (Column 3) and energies of the upper 
level above ground (Colu mn 4). The critical dens ities of each transition 
(Column 5) are taken from lRohlfs & WilsoiJ l2OO0l) . 



Molecule 


AJ 


Vtrans/ GHz 


EJK 


"crit/cm 3 


13 CO 


3-2 


330.58 


31.7 


2 X 10 4 


c 18 o 


3-2 


329.33 


31.6 


2 X 10 4 


H' 3 CO + 


4-3 


346.99 


41.6 


3 X 10 7 



2 OBSERVATIONS 

2.1 Description of observations 

The observations were made over a total of 12 nights between 
6 September and 2 October 2010, using HARP at the JCMT 
(Buck le et al. 1 l2009h . In total, 11.5 hours of data were taken, 
with an area of 400" X 600" (centred on postion 05 h 46 m 40\ 
+00°00'00"(J2000)) being mapped. 

All the standard telescope observing procedures were fol- 
lowed for each night of observations, with regular pointings and 
focusing of the JCMT's secondary mirror. Standard spectra were 
also taken towards various well-known calibration sources, to ver- 
ify that the intensity in the tracking receptor matched recorded 
standards to within calibration tolerance before subsequent obser- 
vations were allowed to continue. The fully-sampled maps were 
taken in raster position-switched mode, an 'on-fhe-fiy' data collec- 
tion method where the HARP array continuously scans in a direc- 
tion paralle l to the sides of the map to produce a full-sampled map 
of the area ( Buckle et al ] |2009h . The telescope is pointed at an off- 
position (05 h 43 m 44 s , +00°21'42.2"(J2000 coordinates)) after every 
row in the map to obtain background values that are subtracted from 
the raw data. A low-resolution 12 CO 7 = 3-2 'stare' observation 
was first performed towards this position to verify that it was indeed 
emission-free. Data are presented in units of corrected antenna tem- 
perature T^, which is related to the main b eam temperature (r mb ) 
using T m b = T^/T] m \, (Rohlfs & Wilson 2000). A value of ^ m b = 0.61 
was used, following Buckle et alj ( 120091) . 

Emission from 3 molecules — 13 CO, C 18 and H 13 CO + — 
was observed and details of the transitions and frequencies are 
given in Table Q] The frequencies and energies quoted are taken 
from LAMBDA, the Leiden Atomic and Molecular Database Q. 
The 13 CO/C ls O data were observed simultaneously, utilising the 
ACSIS high-resolution dual sub-band mode. Each sub-band was 
divided into 4096 61 kHz channels, giving a velocity resolution of 
0.055 km s -1 . The H I3 CO + data was observed in single sub-band 
mode, splitting the 250 MHz bandwidth into 8192 channels, giving 
a velocity resolution of 0.026 km s -1 . The bandwidth of the ob- 
servations is around 200 km s -1 , and has a range of -90 km s _I to 
120 km s -1 . Overall, the weather was good. The mean system and 
receiver temperatures were 330.5 K and 90 K (at 345 GHz) respec- 
tively, and the opacity at 225 GHz ranged from 0.065 to 0.075. 

2.2 Data reduction 

The data were reduced using the Starlink project software, using 
smurf makecube routines to convert from time-series to spectral 

1 http://home.strw.leidenuniv.nl/ moldata/ 

2 http://starlink.jach.hawaii.edu/starlink 



cubes and kappa routines to mask out poorly-performing detectors 
(chpix), remove linear baselines (mfittrend) and alter the image 
resolution (sqorst). 

The final data were sampled onto a 6" per pixel cube, using 
Gaussian gridding kernels with FWHM 9" for all three molecules. 
Taking into accou nt the JCMT beamwidth (14" at 345 GHz 
jBuckle et al]|2009h ), the effective FWHM for the cubes is 17.7" 
for 13 CO and C ls O, and 16.8" for H 13 CO + . The data were also re- 
binned to a velocity resolution of 0. 1 km s~ 1 . At this resolution, the 
mean l<x rms noise on the maps for the three molecules was 0. 10 K, 
0.08 K and 0.03 K for 13 CO, C 18 and H l3 CO + respectively. 



3 RESULTS 

3.1 Integrated intensity and structure 

Figure [2] shows the integrated intensity f 7*^dv images of the 
three molecules (all integrated between 8.5 and 12.5 km s ). The 
SCUB A 450 u.m emission obtai ned from the CADC Legacy Cata- 
logue jPi Francesco et all2008l) is also included for comparison. 

Of the three molecules, the C 18 data most closely resembles 
the optically thin SCUBA dust emission. The large-scale filamen- 
tary structure in the region that is seen in the dust emission is also 
seen in the C 18 0: there is a clear central ridge running E-W which 
contains several smaller dense condensations, and a second large 
filament running almost perpendicular to the ridge in the west of 
the region. Some smaller filaments can also be seen, running diag- 
onally off the main ridge to the southeast. A cavity structure cen- 
tred on 05 h 46 m 40 s , -00°01' can also be seen in the C ls O, and to 
a lesser extent in the 13 CO. The 13 CO emission (being the most 
abundant) covers most of the region, but demonstrates little fila- 
mentary structure. The brightest emission originates in the centre 
of the region, corresponding to the LBS 10 clump. The H 13 CO + 
emission, by contrast, is concentrated in small, dense clumps, with 
no extended emission. Almost all of the emission lies in the cen- 
tral region, although the brightest emission is seen to the west — a 
feature which is conspicuously absent in the CO isotopologues. 

3.2 Velocity structure and kinematics 

The left-hand panel of Figure [3] shows the line centre velocity (as 
determined by a best-fit Gaussian to the 13 CO spectrum at each spa- 
tial pixel) over the entire region. We can see that there is a large- 
scale velocity gradient in a W-E direction over the region — emis- 
sion in the west generally peaks between 9.0 to 9.5 km s" 1 while 
emission in the east generally peaks between 10.5 to 11.0 km s 
— which could indicate differential rotation across the region as a 
whole. 

Sample 13 CO and C ls O spectra from NGC 2068 can be seen 
in Figure [4] and differ greatly depending on the part of the region 
observed — the spectra in the left and right panels of Figure[4]come 
from the east and west of the region respectively. The 13 CO emis- 
sion from the west (corresponding to the LBS 17 clump) has a more 
complex structure — a main peak between 9.5 to 10.0 km s _I and 
either a lower peak between 10.5 to 11.0 km s" 1 or a broad low 
shoulder that extends out to 12 km s" 1 . Double-peaked structure 
can arise in an optically thick line as a result of self-absorption 
or infall, however in this case, the structure is most likely due to 
multiple velocity components along the line of sight. The optically 
thinner C l8 also shows the same multiple-peaked behaviour (with 
more distinctly separated peaks that still occur at similar values 
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RA (J2000) RA (J2000) 

Figure 2. (a) Integrated Intensity of 13 CO (contours from 3 - 30 K km s~' in 3 K km s _I intervals); (b) C ls O (contours from 1 - 9K km s in 1 K km s _1 
intervals); (c) H 13 CO + (contours from 0.1 - 0.5 K km s _1 in 0.1 K km s _1 intervals). Crosses show positions of SCUBA cores identified by M01. (d) 
SCUBA 450 u.m emission obtained from the CADC archives with M01 core positions, and clump positions for LBS 10 and 17, marked with circles and the 
corresponding clump/core number. 




RA (J2000) RA (J2000) 

Figure 3. Left: A velocity map for I3 CO of the line centre velocity at each spatial pixel. Overlaid are 13 CO integrated intensity contours (integrated between 
8.5 - 12.5 km s" 1 ). Right: SCUBA 450 u.m emission, with 13 CO red and blue line wing contours overlaid. Black circles indicate positions of the M01 SCUBA 
cores, and magenta diamonds show protostellar (Class I) cores as determined by NWT from Spitzer IRAC data. 
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Figure 4. Sample (single-pixel) spectra from different parts of the re- 
gion. 13 CO (black line with higher emission) and C ls O (red line) spectra 
from position 05 h 46 m 24.6\ +00°01'06"in LBS 17 (right) and05 h 46 m 46.2 s , 
+00°01'06"in LBS 10 (left) plotted on the same axes. 
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Figure 5. Sample (single-pixel) spectra from different parts of the region. 
C 18 (black line with higher emission) and H 13 CO + (red line) spectra 
from position 05 h 46 ra 31.3 s , -00°02'31.5"in LBS 17 (right) and position 
05 h 46 ra 47.8 s , +00°00'37.5"in LBS 10 (left) plotted on the same axes. The 
H 13 CO + emission is multiplied by 2 times for right-hand plot, and by 5 
times in the left hand plot. 



to the 13 CO), demonstrating the presence of multiple components. 
The emission in the centre and east of the region (corresponding to 
the LBS 10 clump; seen in the left-hand plot of Figure[4]l is much 
simpler, with both 13 CO and C ls O showing single-peaked spectra, 
generally peaking between 10 to 11 km s" 1 . The 13 CO spectra are 
generally fairly wide, with FWHM values of around 2 km s ; the 
C 18 spectra are narrower with FWHM values of around 1 km s" 1 . 
The H 13 CO + spectra (which come almost exclusively from the LBS 
10 clump) are all single-peaked and narrower, generally peaking 
between 10 to 11 km s -1 , and have FWHM values of about 0.5 
km s _1 , as can be seen in Figure[5] 

Figure [TJ shows that the emission in the west (LBS 17) coin- 
cides with an optically-obscured dust lane, while the emission in 
the centre and east (LBS 10) lies on top of the bright reflection 
nebula. This leads to the theory that the emission in the west orig- 
inates from gas and dust in front of the nebula, with the rest of 
the emission originates from deeper in the cloud, behind the neb- 
ula. This theory is supported by polarimetry analysis carried out by 
iMatthews and Wilsorj J2002h . who in addition, argue that while the 
emission originates from two distinct regions, these two regions are 
most likely not completely spatially separated, implying that the fil- 
ament does not lie in the plane of the sky. Instead, the dust filament 
is thought to be twisted with the western edge lying in the fore- 
ground of the optical reflection nebula while the eastern edge lies 
behind it. This orientation of the filaments helps explain the multi- 
ple peaks in the molecular emission — the lower velocity compo- 
nent (that also generally has higher T* h values and is prevalent in 
the western part of the map) lies in the foreground and the higher 
velocity component lies in the background. 



3.3 High-velocity material 

There is evidence of line-wing emission in the 13 CO data (and to a 
much lesser extent, also in the C 18 data), an example of which can 
be observed in Figure|6] In order to capture as much of the outflow 
as possible without contamination from the central emission peak, 
we have integrated from to 8.5 km s -1 for blue-shifted emission, 
and from 12.5 to 20.0 km s _I for red-shifted emission. This inte- 
grated line-wing data is overlaid on SCUBA data in the right-hand 
panel of Figure[3] to demonstrate correlation with the M01 SCUBA 
cores. There are at least 4 separate cases of spatially separated red 
and blue peaks, which we take to be bipolar outflows. 

The isolated bipolar outflow to the west clearly originates 
from the protostellar core at its centre. This core is identi- 
fied as LBS17-MM9 /OriBsmm51, and has also been studied by 
iGibb & Little] feOOOh who concluded that the outflow was driven 
by a deeply embedded Class YSO. The other outflows in the 
centre of the map are much harder to disentangle and match to 
a particular core due to the high density of identified SCUBA 
cores in the central region. LBS10-MM7/OriBsmm35 and LBS10- 
MM3/OriBsmm38 also appear to be at the centre of bipolar out- 
flows, have both been identified by NWT as Class I YSOs, and 
are therefore good outflow source candidates. It is not possible to 
identify them positively as the driving sources of the outflows, due 
to the proximity to LBS10-MM5/OriBsmm3 4 (also identified as a 
Class I YSO) and LBS10-MM6/OriBsmm36. lMitchell et aflfcOOll) 
mapped the region using 12 CO 7 = 3-2 data, but were also unable 
to disentangle the outflows and came to similar conclusions. It is 
likely that since the cores are so close together, without very high 
resolution observations (using an interferometer like ALMA for ex- 
ample), it is not possible to pinpoint the exact core driving each 
outflow. There are patches of both red- and blue-shifted emission 
in the east that are less well-collimated than the others in the region, 
that could originate from LBS10-MM9/OriBsmm39 (identified as 
a Class I YSO). There is also a patch of red-shifted emission in 
the west that could orginate from LBS17-MM6/OriBsmm47, how- 
ever there is no corresponding blue-shifted emission, so we cannot 
postively identify that core as an outflow driving source. 



4 CLUMP DECOMPOSITION 

As we are interested in probing the core kinematics and struc- 
ture, we decompose our C ls O and H 13 CO + emission into smaller 
clumps for analysis, as well as for comparison with numerical sim- 
ulations carried out by others. Both these molecules trace optically 
thin gas at relatively high densities (n ail >10 4 cm~ 3 ), and are there- 
fore likely to reveal information about the dense cores and their 
envelopes/disks. The shapes of the gas clumps produced are also 
less likely to be affected by high-velocity outflows or background 
diffuse gas, as would be the case if we were to perform clump de- 
composition on the 13 CO emission. 

4.1 Data processing and initial results 

We have used the Starlink cu pid implementation of th e gauss- 
clumps algorithm, developed bv lStutzki & Guested <fl99oT) . to iden- 
tify Gaussian clumps in our molecular emission. The algorithm fits 
a triaxial Gaussian profile to the highest peak in the map, then sub- 
tracts the fitted profile from the map, iterating until a user-defined 
termination criteria is reached. We used this method instead of the 
more common clfind algorithm ( I Williams, de Geus. & Blitzl 19941) 
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Figure 6. 13 CO spectra from 3 single pixels around the position of LBS17-MM9 (a known Class YSO), demonstrating the presence of a strong outflow 
source. The red and blue spectra (to the left and right) show the outflow wings clearly. The central pixel still shows residual wings as the outflow is very strong. 



Table 2. Results of the gaussclumps clumpfinding on C ls O and H l3 CO + . 
Column 2: the specified RMS value for each data set, which determines 
the gaussclumps termination criteria; Column 3: the total number of clumps 
found; Column 4: the number of M01 SCUBA cores with a clump match 
(there were 31 M01 cores in total); and Column 5: the number of clumps 
without a core match. 



Molecule 


RMS/ K 


Clumps 


M01 matches 


Noise clumps 


H 13 CO+ 


0.024 


8 


8 (26%) 


1 (13%) 


c 18 o 


0.9 


26 


19 (65%) 


11 (42%) 



as gaussclumps allows more than one clump to be fitted to a partic- 
ular peak in the data, and therefore permits clumps to overlap. This 
is especially useful in the high density centre of the region, as it is 
not always clear where emission from one core ends and another 
begins, and the algorithm used must be good at handling blended 
clumps. 

The termination criteria for the gaussclumps algorithm were 
set to be the following: if more than 10 clumps were fitted with 
peak values less than twice the RMS value; or if 10 consecutive 
failures occurred in the fitting process. To prevent this, the noisy 
edges of the cubes had to be removed and the image cubes were 
cropped in kappa to just encompass the main filaments of emis- 
sions. Any large-scale (> 2.5') background fluctuations in emission 
were removed using cupid findback to smooth the cubes spatially. 

The gaussclumps fits (i.e. the sum of all the clumps found for 
each molecule) are plotted over the integrated intensity observa- 
tions for both molecules in Figure [7] We can see that the gauss- 
clumps algorithm produces clumps that trace the emission fairly 
closely (especially in the case of the H 13 CO + ), with the fitted cubes 
containing around 60% of the original flux. The remaining flux 
is mainly contained in clumps that are smaller than the effective 
beamwidth, and are hence discounted in the final analysis. 

gaussclumps found 8 H 13 CO + clumps and 26 C I8 clumps, 
using the RMS values shown in Table [2] The ellipses shown in 
Figure [8] represent the Gaussian FWHM in the two spatial dimen- 
sions, as calculated by the gaussclumps algorithm. Each clump has 
been matched by eye with the closest M01 core(s) — all cores 
that lie within the ellipse of a clump are matched to that particu- 
lar clump. All but one of the clump-core matches for H 13 CO + pro- 
duced a one-to-one correspondence, indicating that H 13 CO + traces 
the denser core material; the C ls O clumps on the other hand tend 
to be matched to more than one SCUBA core, as they trace a larger 
amount of the less dense material surrounding the cores. 
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Figure 7. gaussclumps fits are overlaid on integrated intensity C O (top) 
andH 13 CO + (bottom) maps of the observations. C ls O contours are between 
1.0 - 4.0 K km s" 1 , at 0.5 K km s" 1 intervals, and then between 4.0 - 10.0 
K km s , at 1.0 K km s _I intervals; H 13 CO + contours are between 0.05 - 
0.5 Kkms" 1 , at 0.05 Kkm s~' intervals. 



4.2 Comparison with simulations 

Smit h et al. IHooi), hereafter SCB09, developed a gravitational po- 
tential clumpfinding algorithm and applied it to simulated data, 
identifying two different types of p-cores (potential cores): bound 
p-cores which quickly collapse and form stars that accrete core 
mass, and composite p-cores which were bound at one point dur- 
ing the simulation but do not necessarily remain bound by the end 
of the simulation. These p-cores are generally smaller than most 
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Figure 8. Top: C ls O integrated intensity, Bottom: H 13 CO + integrated inten- 
sity. Ellipses show the FWHM clump extents as determined by gaussclumps 
for each molecule, and crosses denote the clump centres. Circles show the 
positions of the M01 SCUBA cores. 



dense cores that can be identified observationally, and can only be 
identified when positions and velocities are known in all three di- 
mensions. This is impossible to replicate observationally, and their 
closest observational counterparts are probably high-re solution ob- 
serva tions of cores in very nearby molecular clouds ISmith et al.l 
120091) . Nevertheless, we have found it interesting and useful to 
compare the clumps we identify using more traditional clumpfind- 
ing methods, with those identified using the new gravitational po- 
tential method. 
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R t/ /AU M lir /M f , 

Figure 9. Left: Effective radii for the H 13 CO + and C 18 clumps; Right: 
Calculated virial masses for the H 13 CO + and C ls O clumps. The C ls O 
clumps have a larger range of radii and masses than the H I3 CO + clumps. 

5.1 Shapes and Sizes 

The FWHM sizes (in both the spatial dimensions) for each clump 
were calculated for the two molecules by the gaussclumps algo- 
rithm, and the values were then deconvolved using the effective 
beamwidth as calculated in Section 2.2. The mean deconvolved as- 
pect ratios of the clumps for both molecules are similar: the value 
for H 13 CO + clumps is 0.80 ± 0.15, while that for the C ls O clumps 
is 0.77 ±0.16, where the errors given are the standard deviation^. 
The clumps all tend to follow the orientation of the filament they 
belong to, indicating that they retain the geometry of the filament in 
which they were formed. As shown in Figure [8] for both C ls O and 
H 13 CO + , the major axis is aligned with the filamentary structure 
seen in the integrated intensity map. 

The clump effective radii r cff have been calculated by tak- 
ing the geometric mean of the deconvolved FWHM sizes in the 
two spatial dimensions. The H 13 CO + clumps are on average much 
smaller than the C I8 clumps, as can be seen in the left-hand 
plot of Figure [9] The effective radii of the H l3 CO + clumps is 
3600 ± 900 AU, compared with that of C ls O which have effec- 
tive radii of 6200 ± 2000 AU. The M01 SCUBA cores have a de- 
convolved radii of ~ 2500AU l lMotte et ai]|200ll) . which implies 
roughly a one-to-one correspondence with the H 13 CO + clumps, 
while each C ls O clump will encompass more than one M01 core. 
Comparing our clumps with the bound and composite p-cores of 
SCB09 (which have effective radii of 2400 AU and 3700 AU re- 
spectively), we find that the H I3 CO + clumps have r eff comparable 
to the average p-core radius, and are therefore more likely to corre- 
spond to individual cores. The C 18 clumps on the other hand are 
much bigger and encompass multiple p-cores; these are expected to 
correspond to the bound p-cores which have masses that are about 
half that of the C 18 clumps on average (see Table|4] 



5 CLUMP PROPERTIES 

The C ls O is expected to trace gas at lower densities than the 
H 13 CO + , and the two molecules will therefore probe different re- 
gions of star-forming cores. In order to investigate these differ- 
ences, we examine and comp are the properties o f the clumps pro- 
duced for the two molecules. Ilkeda et alj 12009J) also conducted a 
survey of dense gas cores in Orion B (which included NGC 2068), 
using H 13 CO + 7 = 1-0 molecular emission, and identified 151 
dense cores using the clfind algorithm l lWilliams. de Geus. & Blitj 
1 19941) . We present a comparison of our clump properties with 
both SCB09 core properties, and with those cores identified by 
Ilkeda etalj d200St) in Table 4. 



5.2 Linewidths 

The average linewidth for each clump calculated by gaussclumps — 
the FWHM of the Gaussian fit to the average spectrum for the 
clump — is generally smaller than the linewidth at the clump cen- 
tral position. In fact, this discrepancy is so much greater for the 
C 18 clumps than for the H 13 CO + clumps that the clump average 
linewidths for C 18 are smaller than those of H I3 CO + : 0.18 ± 0.05 
km s" 1 for C 18 compared to 0.24 ± 0.08 km s" 1 for H 13 CO + . 
However, when a Gaussian fit to the spectra at the clump centre 
is performed, an average value of 0.28 ± 0.02 km s" 1 for C 18 

3 In all cases, unless otherwise stated, all errors quoted are the standard 
deviations of the sample 



© 2012 RAS, MNRAS OOO.mTRl 



8 S. L. Walker-Smith et. al. 



I 1 H^CO* 

c:: c"o 



0.10 0.15 20 25 30 0.35 40 0.45 0.50 
Linewidths/kms 1 




Internal ID Velocity Dispers 



Figure 10. Left: Histogram showing the linewidths of the clumps in the 
C ls O (dashed line) and H 13 CO + (solid line) datasets. Right: Histogram 
showing the ID internal velocity dispersions of the p-cores in the bound 
(solid line) and composite (dotted line) datasets, taken from SCB09. 



and 0.27 ± 0.04 km s" 1 for H 13 CO + is obtained. These latter val- 
ues of C I8 are a better reflection of the actual observed spectra 
of the molecules (see Figure [5}, where the C l8 emission line is 
visibly wider than that of H 13 CO + . This result appears to contra- 
dict the well-kno wn result that the linewidth increases with the size 
of region traced iLarsoiJll98lh . However, this can be explained by 
the fitting mechanism in gaussclumps: the central peaks are gener- 
ally well-defined and easy to fit; but at the clump edges, the emis- 
sion decreases and the fits are made to peaks that have much lower 
SNR and tend to be much narrower. Therefore the average spectra 
produced have a significant contribution from very narrow spectra, 
causing them to be narrower than the peak spectrum. This effect is 
more noticeable for the C ls O clumps than for the H I3 CO + clumps, 
because the latter tend to have a much sharper cutoff as the emis- 
sion is very concentrated. Therefore, in all subsequent calculations, 
to standardise the comparisons and to avoid biasing the spectra due 
to noise contributions, we will use the clump centre linewidths o"id 
and central velocities v c , rather than the average values produced 

by GAUSSCLUMPS. 

The linewidths are shown in Figure [10] with the internal ID 
velocity dispersion of the p-cores from SCB09 for comparison. We 
see that the distribution of both sets of clumps peaks at a much 
higher value than the SCB09 p-cores — SCB09 calculate cr 1D val- 
ues of 0.16 km s _1 and 0.23 km s -1 for bound and composite p- 
cores respectively. This is to be expected for our C 18 clumps as 
they are much larger than the SCB09 cores, and trace more of the 
turbulent diffuse gas surrounding star-forming cores. The H 13 CO + 
clumps trace very dense gas, and would be expected to have much 
smaller linewidths. However, over half of the clumps are associ- 
ated with centrally condensed protostellar objects; these clumps 
exhibit higher linewidths than those that are starless and there- 
fore increase the average clump linewidth. The average linewidth 
of those H 13 CO + clumps not associated with a protostar is 0.17 ± 
0.05 km s _1 , which matches the values obtained for the SCB09 
cores much better. It should be noted that the thermal width of 
H 13 CO + at 16 K is 0.07 km s -1 , making the thermal contribution 
to the linewidth less than 10%. 

We analyse the clump linewidths in greater detail, and also 
discuss our results in relation to theories of star formation, as well 
as timescales and size scales of turbulence, in Section 5.6. 



5.3 Virial Masses 

The virial mass of a spherically symmetric clump 
can be calculated using the f ollowing equation 
dMacLaren. Richardson. & Wolfendaielll988l) : 



Table 3. Various properties of H CO* and C O used in this paper to cal- 
culate column densities and gas masses. Values of Tq, the energy of the 
lowest state, in K, ar e obtained from the Leiden Atomic and Molecular 
Database (LAMBA) ISchoier et alj|2005h. The r e lative abun dances X gas 
were obtained from iFrerking. Lange r. & Wilsorj i 19821) and llkeda et al] 
<2009l) . 



Molecule 


r /K 


Xgas 


H 13 CO + 


4.16 


4.3 x 10"" 


c 18 o 


5.27 


1.7 x 1(T 7 



h r cS crl 



(1) 



where a\ D 



3fX 2 D , 



T e ff is the effective radius of the clump (as 
discussed in Section 5.1), and G is the gravitational constant; k\ is 
a constant whose exact value depends on the density distribution 
of the clump as a function of distance from the clump centre. For 
simplicity we assume that the clumps follow p oc r" 1 density distri- 
butions, yielding k\ = 1. 

The virial mass oftheH 13 CO + clumps varies beween 0.3 - 2.0 
M , while that of the C 18 clumps has a range of 0.2 - 2.7 M (see 
Figure[9). In general, the average C ls O virial mass of 0.86 ± 0.13 
M G is similar to that of the H 13 CO + clumps (0.97 ± 0.21 M Q ). 



5.4 Gas Masses 

The gas mass of the clumps can be calculated from the column den- 
sity of the gas, assuming Loca l Thermal Equilibrium (LTE). Fol- 
lowing |RohS^wiioJl|2003) and assuming values for the differ- 
ent molecules as presented in Table [3] 



N r 



5.70 X 10' 



8.70 x 10 I6 nT 



exp(-41.7A77; 
T 

1 I'X 



(2) 



(3) 



exp(-31.7tf/7; x ) 

The output from cupid — the 'clump sum' — is the sum of the 
data enclosed within the gaussclumps ellipse (defined by the Gaus- 
sian FWHM). The integrated intensity J T* A dv is obtained by mul- 
tiplying this clump sum by the velocity resolution of the cubes (0. 1 
km s" 1 ). In both cases, we take the excitation temperature to be 
16.1 ± 4.7 K, the average value calculated for the entire region by 
llkeda et al.l | |2009|) . using NH 3 transitions. Using a single tempera- 
ture to represent all the different cores will introduce a certain in- 
accuracy to the gas masses calculated, as it will not necessarily re- 
flect the different conditions within each clump. However, it is the 
least biased method we could choose and allows consistent com- 
parisons for all clumps identified. Additionally, our most optically 
thick molecule 13 CO is only barely optically thick (with opacity 
values of t ~ 0.9), and will not give accurate excitation tempera- 
tures, tending towards underestimation; we have therefore chosen 
not to use excitation temperatures calculated from 13 CO. 

We obtain the masses of the clumps by substituting represen- 
tative values for the two molecules into the following equation: 

M gas = D 2 ( Ao- A/?)/jh2 m H2 X ga ' s N gas (4) 

We take the distance to the cloud D = 415 pc and adopt a mean 
molecular weight per H2 molecule of /jh2 = 2.72 to include helium; 
we also use a pixel size (Ao- and A/3) of 6" and take the abundances 
of H 13 CO + and C ls O relative to hydrogen (X H n CQt and X c k ) to 
be as shown in Table[3] 
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Figure 11. Plots of (Top) M vn . against M gas and (Bottom) M v ; r against 
M imt , for both H 13 CO + (red) and C 18 (blue). The M vir = 2M gas (bound) 
and M v ; r = M gas (equipartition) lines are also indicated using black and 
grey lines respectively. 

The values of the gas masses calculated using this method 
range from 0.5 to 1.4 M Q for H 13 CO + , and from 0.6 to 12 M Q for 
C 18 0. As expected, the C I8 clumps contain more gas mass due to 
their much larger sizes, having an average mass of 3.7 ± 3.0 M Q 
compared to the average mass of 0.92 ± 0.24 M G calculated for 
H 13 CO + . 

5.5 Bound Ratio 

One method of separating and identifying clumps that are prestellar 
and likely to collapse to form protostars rather than transient struc- 
tures, is to determine whether they are bound. A fairly straightfor- 
ward method of doing this is by comparing their virial masses (a 
measure of their internal energy) with their gas or dust masses (a 
measure of the potential energy in the clump). Clumps are generally 
considered to be bound if M vir < 2M gas ; however, the lower limit 
for equipartition, where the clumps are confine d but not necessar- 
ily se lf-gravitating, is taken to be M v i r = M gas l lBertoldi & McKeel 
[1991) . 

Figure QT] (top) shows a plot of clump virial mass against gas 
mass, where all of the C ls O clumps and a third of the H 13 CO + 



clumps lie to the right of the 'bound' line, and all but one of 
the H 13 CO + clumps are at least in equipartition. This implies that 
the larger C I8 clumps are more bound than the smaller H 13 CO + 
clumps. H 13 CO + traces very dense gas and we would expect the 
structures identified to be more tightly bound than the lower density 
structures traced by C ls O, making our result appear contradictory. 
For comparison, we also plot clump virial masses against their dust 
masses (see Figure|TT] bottom), which are calculated from the dust 
masses of the corresponding M01 SCUBA cores. If (as is usually 
the case for C ls O) more than one M01 core corresponds to a partic- 
ular clump, then the core masses are summed. (Those gas clumps 
without an M01 core match are omitted in this part of the analysis.) 
It is reassuring to see that the majority of the H I3 CO + clumps now 
lie to the right of the 'bound' line, and that the only clumps lying to 
the left of the equipartition line (and therefore unbound) are C 18 
clumps. 

The dust masses seem to give a more accurate reflection of 
the mass within the H 13 CO + clumps probably because they are of 
similar sizes, while giving an underestimation of the mass in the 
C ls O clumps, which are much bigger and encompass more of the 
envelope material surrounding the dense cores. However, there are 
sources of error (for example, spatial filtering in dual-beam switch 
observing mode) that could bias the dust masses towards lower val- 
ues. In addition, there are intrinsic sources of error and uncertainty 
present, such as dust emissivity and the dust/gas mass ratio, as well 
as contamination from molecular gas like 12 CO, that could affect 
dust mass estimates. 

One of the reasons for the C 18 clumps appearing more bound 
is an overestimation of the gas mass. C ls O has a lower critical den- 
sity than H 13 CO + and therefore will exhibit a contribution from 
the less dense surrounding envelopes of pre- and protostellar cores. 
This would account for the much larger sizes and higher gas masses 
of C 18 over H 13 CO + — the H 13 CO + clumps represent the denser 
core regions while the C 18 clumps represent both core envelopes 
and centres. In addition, the H 13 CO + may have been sub-thermally 
excited which would result in an underestimation of the column 
density and hence the mass present in the clump — this theory is 
supported by the fact that a large proportion of the H 13 CO + clumps 
have higher dust masses than gas masses. 

There are also many sources of error present in the calcula- 
tion of gas masses, the most significant of which is the value of 
the abundance ratios of H 13 CO + and C 18 0, which can cause the 
calculated gas masses to vary by a factor of 3. 

Finally, the comparison of virial and gas/dust masses merely 
tells us if the material in the clump is bound with respect to itself, 
but does not tell us if the clump is bound with respect to the sur- 
rounding gas environment. There are other (less easy-to-quantify) 
methods of confinement for clumps (e.g. magnetic fields, external 
pressure) that could mean that a clump is in fact a bound (and hence 
non-transient) structure. In fact, converting a linewidth straight to 
a virial mass without accounting for surface pressure will overesti- 
mate the enclosed mass of the core and h ence make it app ear less 
bound than it actually is (for example, see lDib et al.l d2007l) ). 

5.5.7 Overall region boundness 

Although the comparison of virial mass with gas mass may not 
be the most suitable method for determining if small sub-structures 
within molecular cloud regions are bound, it should still give a good 
idea of the condition of the region as a whole. The 13 CO emission 
lines are resolved across the whole region, and will therefore give 
a good estimate of the virial mass of the region; the C 18 emission 
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Table 4. Comparison of results from this paper with referenced literature values. The typical clump mass values used in the case of this work, are the gas 
masses of the clumps, calculated from the column density values. The errors on the last digit are given in brackets for the <xirj values. 



Author 


Region 


Clump type 


Typical Clump Size/ AU 


Typical Clump Mass/M Q 


Mean <xirj/km s 1 


Ikeda et. al (2009) 


Orion A 


H 13 CO + (J=l-0) 


2.9 X 10 4 


2.1-81 


0.52 


Smith et. al.(2009) 


simulated 


p-cores (bound) 


2.4 X 10 3 


0.2-5 


0.16 


Smith et. al.(2009) 


simulated 


p-cores (composite) 


3.7 X 10 3 


0.2 - 10 


0.23 


This work 


Orion B 


H 13 CO+ (J=4-3) 


3.6 X 10 3 


0.5-1.4 


0.27(4) 


This work 


Orion B 


C ls O (J=3-2) 


6.2 X 10 3 


0.6-12 


0.28(2) 



also covers the majority of the region, is optically thin throughout 
(unlike the 13 CO), and should give a good estimate of the amount 
of gas contained within the region. We can thus calculate the virial 
and gas mass of the NGC 2068 region, and determine whether it is 
bound as a whole, using a similar method to that described previ- 
ously. Assuming that the region is spherically symmetric, we there- 
fore calculate the virial mass from: 

M vir =^^k (5) 

The value of r cond — the radius of the NGC 2068 region — is taken 
to be 4' (which translates to ~0.5 pc at a distance of 415 pc); and 
o" tot = 0.61 km s _I , as obtained from the dispersion of 13 CO line 
centre velocities over the region. Substituting these values, M vil . is 
approximately 130 M G . 

Using Equations [3] and [4] we calculate M gas of approximately 
230 M Q from the C' s O emission, which is almost twice the calcu- 
lated virial mass. This indicates that the cloud is most likely in the 
process of global collapse, unless additional support (e.g. magnetic 
fields) were present. In order to provide this support, the magnetic 
field strength would have to be ~ 59^G (follow ing virial equilib- 
rium calculations by Bertoldi & McKe^ J 1 992h "). which is similar 
to the value of ~ 36pG they have measured in Orion B. We should 
not however place too much weight on this value due to the com- 
plex dynamics and irregular shape of the NGC 2068 cloud. 

5.6 Velocity dispersions 

The kinematics of fragmenting cores can be used to investi- 
gate the initial conditions of star formation. An examination of 
the internal velocity dispersions of these cores tells us the mag- 
nitude of turbulence present w ithin the cores on small scales 
( iKirk. Johnstone. & Tafallal2007l) . The relative movement between 
cores, as well as between the cores and their surrounding gas, could 
shed light on the timescale of turbulent decay present in the region 
— whether the turbulence decays on dynamical timescale s, or is 
driven and lasts for longer timescales. lOffner et al .1 J2008h found 
that comparing the core-to-core velocity dispersions of protostel- 
lar and starless cores allowed them to differentiate between their 
driven and decaying turbulence simulations. 

J. 6. 1 Clump internal velocity dispersions 

The o"i D values, shown in Figure [T0l are equivalent to the clump 
internal velocity dispersions and shed light on the turbulent nature 
of the clump interiors. 

We compare the cr 1D values for starless and protostellar 
clumps in Table [5] to determine if there is a statistically significant 
difference between the two types of objects. We define a particular 



M01 core as being protostellar if it has an IRAC Class I match in 
NWT, or if M01 have designated them as protostellar. Our clumps 
are then deemed to be protostellar if a protostellar M01 core lies 
within the extent of the ellipse defined by the gaussclumps algo- 
rithm. We also narrow our range of starless clumps down to only 
include those that are associated with an M01 core, all others are 
deemed to be 'noise' or transient clumps. Using these definitions, 
we find 4 protostellar and 3 starless H 13 CO + clumps; and 5 pro- 
tostellar and 10 starless C 18 clumps. We find that in the case of 
H 13 CO + , the 0"id values are much lower for starless clumps than 
for protostellar clumps; the C ls O cr i0 values also show the same 
trend, but to a lesser extent. 

We performed a Kolmogorov-Smirnoff (KS) test, to determine 
if the two protostellar and starless datasets were drawn from the 
same underlying distribution. The KS statistic (or p-value) gives 
a value that is 1 minus the confidence level with which the null 
hypothesis (that the two samples originate from the same distribu- 
tion) may be rejected. Generally, if the p-value is < 1%, we can 
say confidently that the two samples originate from different un- 
derlying distributions. KS tests give p-values of 15% and 26% for 
H 13 CO + and C ls O respectively, when comparing protostellar and 
starless clump distributions. Therefore at the 85% confidence level, 
the H 13 CO + protostellar and starless clumps originate from differ- 
ent underlying distributions, justifying our clump classification sys- 
tem. There is a greater difference between the two distributions in 
the case of H 13 CO + as it traces the denser gas in the cores, which 
should be more quiescent and also contain less centrally condensed 
gas for starless cores than protostellar cores; the C ls O on the other 
hand contains more of a contribution from the core envelopes (as 
mentioned previously) making the distributions less distinct. 

The tr 1D values are also compared to the sound speed c s in 
the region, to determine how turbulent the clump interiors are. We 
define three turbulent regimes: subsonic (cr iD < c s ), transonic (c s < 
<Tid < 2c s ) and supersonic (2c s < o"i D ). The value of c s is calculated 
using the following formula: 

c s = (6) 
V m 

r ex = 16.1 K is the value adopted earlier in Section 5.4 , and m = 
2.33w?p, assuming a ratio of 1 Helium for every 5 IU molecules, 
giving a value of c 8 = 0.24 km s -1 . We also calculate the bulk gas 
velocity dispersion <x g = 0.61 ± 0.29 km s (from the l3 CO gas 
over the region), to compare the clumps' internal turbulence with 
their surrounding natal gas environment. 

We find that the <Tid values for all clumps are transonic, except 
for the starless H 13 CO + clumps which are subsonic. This is con- 
sistent with the starless clumps corresponding to local minima of 
turbulence, as they are much less turbulent internally than the sur- 
rounding gas. In addition, we see that the <tid values for all clumps 
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Table 5. Mean internal ID velocity dispersions of the Protostellar (P) and Starless (S) C ls O and H 13 CO + clumps; as well as for the two molecules without 
differentiating clump type. The errors on the last digit (standard deviations on the sample) are given in brackets after each value. 



Dispersion 


H I3 CO+ (P) 


C I8 0(P) 


H 13 CO + (S) 


C ls O(S) 


H I3 CO + (All) 


C 18 (All) 


trio/km s -1 


0.30(8) 


0.35(6) 


0.17(5) 


0.30(5) 


0.27(10) 


0.28(8) 


0"id/c s 


1.4(4) 


1.47(23) 


0.86(22) 


1.25(22) 


1.10(4) 


1.2(3) 




0.55(17) 


0.58(10) 


0.34(9) 


0.49(9) 


0.45(17) 


0.46(2) 



are much smaller than the bulk gas velocity dispersion. These re- 
sults are consistent with the predictions of gravoturbulent fragmen- 
tation — i.e. that the protostellar cores are found at the intersection 
points of turbulent shocks and are therefore much less turbulent in- 
ternally than their surrounding natal gas. However, we do note that 
most star formation scenarios agree that the non-thermal support 
should decrease in star-forming cores, so this is not conclusive ob- 
servational evidence of gravoturbulent fragmentation. 



5.6.2 Relative core- gas motions 

The interclump velocity dispersion cr v that we calculate is the dis- 
persion (or standard deviation) of the line centre velocities v c at 
the clump centres, for all the clumps in the region. Initially, with- 
out distinguishing between protostellar and starless clumps, we find 
that the interclump dispersions decrease as the density of gas traced 
by the molecule increases (see Figure [72]): 0.39 ± 0.05 km s _1 and 
0.28 ± 0.08 km s~' for C I8 and H 13 CO + respectively. A KS test 
comparing <x v for the two molecules indicates with 85% confidence 
that they are traci ng gas from different underlying distributions. 
I Andre et al. I J2007t) calculated a core-to-core dispersion of ~ 0.25 
km s" 1 for protostellar cores in the L1688 region of Ophiuchus us- 
ing N2H + , which traces similar densities of gas and is very close 
to the value we calculate for H 13 CO + clumps, but only about two- 
thirds of that for the C 18 clumps. This discrepancy could be ex- 
plained however by the difference in global properties of the overall 
molecular cloud. Orion is a much more turbulent, clustered star- 
forming environment than Ophiuchus and the differences between 
bulk gas motion and individual star-forming cores might be more 
pronounced. The C 18 clumps cover much larger areas and would 
be more likely to be influenced by the bulk gas movements, whereas 
the H 13 CO + clumps are tracing denser regions and are more likely 
to reflect the dispersions between star-forming cores. 

Using the same selection criteria as before, we compare the in- 
terclump dispersions between protostellar and starless clumps. The 
interclump dispersions were also compared with the region sound 
speed and the region bulk-gas velocity, in the same manner as in 
the previous section. The results are presented in Table[6] 

All clumps show transonic behaviour (except the starless 
H 13 CO + clumps which are subsonic), indicating that there is lit- 
tle discernable movement between clumps. Comparing the inter- 
clump dispersions with the bulk-gas velocity dispersion, we find 
that the motion of the clumps compared to their surrounding gas 
is fairly low, with o~ v /o~ g values of 0.46 and 0.80 for H 13 CO + and 
C 18 respectively. This means that the clumps are in some way still 
coupled to the surrounding gas from which they formed, and indi- 
cates that the clumps should be able to move together with their 
surrounding gas under the same gravitational force s, one of the key 
requi rements for competitive accretion to occur dBonnell & Bate! 
l2006h . The relative clump-gas velocity for H 13 CO + is lower than 
that for C ls O, which is most likely due to the C ls O tracing more of 



the surrounding envelope that is expected to be more turbulent, and 
henc e will retain a g reater degree of the movement of the natal gas. 

lOffner et alJ d2008l) compared core-to-core velocity 
dispersions for observational data taken in the Perseus 
Kirk, Johnstone, & Tafallal 120071) and p Ophiuchus regions 



Andre et alJ[2007l) . with those produced in simulations of driven 



and decaying turbulence. They found that the driven run was 
in better agreement with Perseus while the decaying run gave a 
better match for Ophiuchus. They did however note that the two 
simulations were statistically very similar, and that in general the 
distribution of core central velocities did not appear to depend on 
the details of the turbulence. When they differentiated between 
prestellar and protostellar cores in their simulations, they were able 
to distinguish between the two turbulent scenarios. They found 
that in driven turbulence simulations, the prestellar core-to-core 
velocity dispersions were larger than those for protostellar cores; 
conversely for decaying turbulence simulations, the protostellar 
core-to-core velocity dispersions are larger than those for the 
prestellar cores. In both cases, the difference is about a factor of 
1.3. 

Comparing interclump dispersions for our protostellar and 
starless clumps (from the same molecule), we find that H 13 CO + 
protostellar clumps have a higher <x„ than the starless clumps, in- 
dicative of decaying turbulence; the opposite is true for C 18 0, 
where <x„ is higher for starless clumps than protostellar clumps. 
From KS tests, we cannot rule out the possibility that the pro- 
tostellar and starless clumps originate from the same underlying 
distribution (p-values of 26% and 54% for H 13 CO + and C I8 re- 
spectively). The C 18 gas may also be more affected by other fac- 
tors such as multiple velocity components along the line-of-sight or 
outflows, which could all potentially confuse the line centre veloc- 
ity. We therefore believe that the H 13 CO + will give a better picture 
of the star-forming cores themselves, given that they are tracing 
the denser gas. In comparison, c ore-to-core velocity dispersions o f 
N2FP observations in Perseus ( Kirk . Johnstone. & T afalla 120071) . 
provide support for driven turbulence — they calculated a higher 
o~ v for their prestellar cores than for their protostellar cores. As their 
emission traces similar densities to the H 13 CO + in this analysis, the 
type of turbulence present in Perseus is therefore different from that 
present in NGC 2068. However, as we only have a small sample of 
H 13 CO + clumps, we cannot say for certain whether decaying turbu- 
lence is dominant in our region; a larger statistical sample (prefer- 
ably over different regions, to confirm environmental effects) would 
be required to confirm the results. 



5.6.3 Comparison with Decaying Turbulence Numerical 
Simulation 

The main numerical simulation that we are comparing our data 
with, is that of SCB09. The simulation consists of a 3-dimensional 
cylinder (3 x 3 x 10 pc) containing 10 4 M Q concentrated at one 
end, so the top is over-bound and the bottom is under-bound. The 
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Table 6. Interclump velocity dispersions of the Protostellar (P) and Starless (S) C ls O and H 13 CO + clumps; as well as for the two molecules without differen- 
tiating clump type. The errors on the last digit(s) are given in brackets after each value, and are calculated using the formula from lAndre et alj 120071) . 



Dispersion 


H 13 CO + (P) 


C 18 0(P) 


H 13 CO + (S) 


C ls O(S) 


H 13 CO + (All) 


C ls O (All) 


<r,,/km s 


0.27(11) 


0.35(12) 


0.22(11) 


0.49(11) 


0.28(8) 


0.39(5) 


o-vlc s 


1.14(3) 


1.45(3) 


0.92(3) 


2.03(3) 


1.17(2) 


1.61(2) 


<r v /a g 


0.45(14) 


0.57(18) 


0.36(11) 


0.80(25) 


0.46(14) 


0.63(19) 



turbulence in this simulation is not driven but allowed to decay once 
self-gravity is applied; however it is replenished by the release of 
kinetic energy as the molecular cloud collapses. Further details of 
the simulation can be found in the SCB09 paper. A catalogue of p- 
cores (found via a gravitational potential clumpfinding algorithm, 
as mentioned in Section 4.3) is produced at snapshot intervals of 0. 1 
fdyn (4.7 x 10 4 years). For each p-core, the catalogue contains in- 
formation about the centre-of-mass (CoM) velocities and positions, 
and their masses, as well as their potential, kinetic and thermal en- 
ergies. In their simulation, SCB09 have defined a subset of 'bound' 
p-cores, for which E mt > 1, where E mi is defined as: 



I I Region A 
CZ1 ReglonB 



l|£ P 



+ E K 



(7) 



and E lherm is the thermal energy of the clump, E K is the kinetic en- 
ergy calculated with respect to the CoM velocity of the clump and 
E p is the potential energy of the clump calculated from the depth of 
the background- subtracted potential well. These p-cores are there- 
fore bound with respect to the environment in which they were 
formed, and not merely in isolation (as is the case for the clumps 
we identified using traditional clumpfinding methods). Bound p- 
cores are assumed to be pre-stellar in this simulation, as they go 
on to form sinks that accumulate mass and serve as the observable 
manifestation of star formation. As such, we would like to make 
velocity dispersion comparisons of our H 13 CO + and C 18 clumps 
with these pre-stellar cores. 

We use p-core catalogues from 4 snapshots in time: 1.025, 
1.050, 1.125 and 1.180 fd yn , corresponding to a span of time be- 
tween 4.7 x 10 5 and 5.6 x 10 5 years. We combined the 4 catalogues 
into a single data-set, to eliminate any time-dependent effects and 
to increase the sample size. A KS test performed on the 4 snapshots 
gave p-values > 20 %, implying that they originate from different 
underlying distributions, and we are justified in assuming the sam- 
ples to be independent when combining them. We also refine our 
sample by only choosing bound p-cores, using the same criterion 
as SCB09, that E rat > 1. Our final dataset consists of 533 bound 
prestellar p-cores, which is large enough to be statistically mean- 
ingful. 

Initially, we calculate the intercore velocity dispersion cr r of 
all the p-cores in the simulation, in the same way as we did be- 
fore. We find that cr v = 0.602 + 0.017 km s , which is much larger 
than the individual interclump dispersions of the H 13 CO + and C ls O 
clumps. 

To better match the size of NGC 2068, we then extract 2 
parsec-sized regions, which are taken from different parts of the 
simulation, and so have differently-bound initial environments. Re- 
gion A correponds to the initially less-bound environment, while 
Region B corresponds to the initially more-bound environment. 
Figure Q/2] shows a comparison between the observed and simu- 
lated intercore velocity dispersions. We calculate cr v for the p-cores 
within each region, and obtain values of 0.41 ± 0.05 km s _1 and 
0.36 ± 0.05 km s" 1 for Re gions A and B respectively. The values 
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Figure 12. Histograms showing the spread of clump line centre velocities 
for the H 13 CO + and C ls O clumps (left); and the equivalent spread of line- 
of-sight velocities calculated for the bound p-cores in SCB09 (right). 



calculated for the smaller regions are very close to the interclump 
dispersions for the C 18 clumps (see Table|6]l, which makes sense 
as the size scales over which we are calculating the dispersions 
are similar. The H 13 CO + clumps are concentrated in a smaller area 
than the C 18 clumps, and this is reflected in their interclump dis- 
persions, which are 70% of the p-core cr v values. We therefore sug- 
gest that when calculating and comparing interclump dispersions, 
the size of the region considered is an important factor, as the <x„ 
appe ars to vary wi th region size. This is similar to the Larson rela- 
tion jLarsonlll98ll) of cr(km s~')= l.lOL(pc) 038 , which reflects the 
underlying turbulence scaling between cr v and region size. 

We find that the simulation — which models a decaying turbu- 
lence environment — is able to reproduce quite closely the effects 
and properties that are seen observationally, when comparing sim- 
ilar size-scales. However, there is still insufficient evidence to say 
conclusively that the region we observe is dominated by dynami- 
cally decaying turbulence. 



6 CONCLUSIONS 

We have presented 13 CO, C 18 and H 13 CO + data for the NGC 2068 
region in Orion B, the combination of which allows us to probe 
different densities and size-scales, and determine the structural and 
kinematic properties of the region. 

The 13 CO is widespread and detected over the entire re- 
gion, while C I8 has a more filamentary structure and follows the 
SCUBA dust emission much more closely; the H 13 CO + is confined 
to small clumps, the majority of which are found in the region cen- 
tre. The 13 CO and C ls O show evidence of multiple velocity com- 
ponents along the line-of-sight, especially in the west of the region; 
the H 13 CO + remains single peaked throughout, with a line centre 
velocity of 10.86 ± 0.10 km s _1 , reflecting the systemic velocity of 
the region. The 13 CO also shows evidence of several bipolar out- 
flows in the region, and we have matched these to SCUBA cores, 
as far as possible given the density of sources in the central region. 

We have used the gaussclumps algorithm to decompose the 
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C 18 and H 13 CO + data into Gaussian clumps, and calculated var- 
ious clump properties, including size/shape, linewidths, and virial 
and gas masses. The C l8 clumps have higher linewidths than the 
H 13 CO + clumps on average, and are significantly larger than the 
H 13 CO + clumps. We therefore conclude that the H 13 CO + traces the 
denser inner parts of the star-forming cores, while the C 18 likely 
traces more of the less dense, outer envelopes of the cores. Using 
the ratio of virial to gas mass, we have found that the C 18 clumps 
mostly appear more bound than the H 13 CO + clumps. However this 
may not be the best method of calculating how bound a particu- 
lar clump is — the C 18 clumps appear more bound as there is a 
greater contribution from the surrounding envelope; using the dust 
mass instead of the gas mass gives an alternative view of the mass 
contained within the H I3 CO + clumps, making them appear much 
more bound. 

We have classified our molecular clumps as protostellar or 
starless based on a match with a known YSO or SCUBA core 
respectively. We have also used the internal and interclump ve- 
locity dispersions for the two molecules to investigate star for- 
mation mechanisms. The internal velocity dispersions of the star- 
less H 13 CO + clumps (probably the most likely objects to repre- 
sent prestellar cores) are subsonic, suggesting that the cores were 
formed at local turbulence minima, consistent with gravoturbulent 
fragmentation. The low values of the interclump velocity disper- 
sions compared to the bulk gas dispersion (as traced by the 13 CO) 
also indicate a certain degree of coupling between the possible star- 
forming cores and their natal gas environment, fulfilling one of 
the key requirements for competitive accretion to occur. The in- 
terclump velocity dispersions match numerical simulations of de- 
caying turbulence, when similar sized regions are compared; the 
H 13 CO + protostellar and starless interclump velocity dispersions 
also match predictions for decaying turbulence well. However, due 
to the small size of our sample, we cannot say for certain that decay- 
ing turbulence is dominant in the region. A larger statistical sample 
(over regions with different gas environments) would be required 
to conclusively determine the type and decay timescale of the tur- 
bulence present in this region, and how the nature of the turbulence 
varies between regions. 
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APPENDIX A: CLUMP CATALOGUES 

A catalogue of the clumps found by gaussclumps for the H 13 CO + 
and C 18 data is provided in this appendix, in Tables [ATI and lA2l 
respectively. 
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Table Al. Catalogue of gaussclumps clumps found for H 13 CO + emission in NGC 2068. Column 2: Deconvolved FWHM size of each clump in the 2 spatial 
dimensions; Column 3: Peak intensity in each clump; Column 4: Internal velocity dispersion crirj at the peak of each clump; Column 5: Line centre velocity 
at the peak of each clump; Column 6: Ratio of the FWHM sizes of each clump; Column 7: Effective radius of the clump; Column 8: Virial mass of the clump; 
Column 9: The sum of values extracted from within the ellipse defined by the clump FWHM; Column 10: Mass of clump calculated using the values in column 
9 assuming LTE; Column 1 1 : Ratio of the gas mass to the virial mass; Column 12: Mass of the M01 SCUBA cores associated with each clump. 



Clump 


Size 


Peak 


C"1D 




Axis 


'oil' 


M vir 


Clump 




Bound 


Msmm 


No. 


/(AU X AU) 


/K 


/km s 


/km s -1 


ratio 


/pc 


/Mq 


sum /K 


/Mo 


ratio 


/Me 


1 


2200 x 2100 


0.55 


0.28 


10.54 


0.99 


0.022 


0.73 


53.6 


0.93 


1.3 


3.00 


2 


2600 x 8200 


0.50 


0.19 


11.05 


0.63 


0.037 


0.42 


63.8 


1.10 


2.6 


1.05 


3 


1 1000 x 3300 


0.42 


0.22 


10.94 


0.55 


0.044 


0.89 


48.1 


0.83 


0.9 


3.50 


4 


3500 x 3100 


0.42 


0.43 


10.93 


0.97 


0.026 


1.63 


53.5 


0.92 


0.6 


1.35 


5 


8900 x 5500 


0.35 


0.16 


11.17 


0.75 


0.046 


0.34 


80.5 


1.39 


4.1 


0.40 


6 


5000 x 2400 


0.26 


0.42 


10.40 


0.82 


0.029 


2.04 


51.0 


0.88 


0.4 


1.35 


7 


5200 x 3200 


0.22 


0.21 


11.14 


0.85 


0.031 


0.80 


40.3 


0.70 


0.9 




8 


7400 x 5400 


0.22 


0.27 


10.74 


0.84 


0.042 


0.93 


33.4 


0.58 


0.6 


0.90 



Table A2. Catalogue of gaussclumps clumps found for C ls O emission in NGC 2068. Column 2: Deconvolved FWHM size of each clump in the 2 spatial 
dimensions; Column 3: Peak intensity in each clump; Column 4: Internal velocity dispersion crirj at the peak of each clump; Column 5: Line centre velocity 
at the peak of each clump; Column 6: Ratio of the FWHM sizes of each clump; Column 7: Effective radius of the clump; Column 8: Virial mass of the clump; 
Column 9: The sum of values extracted from within the ellipse defined by the clump FWHM; Column 10: Mass of clump calculated using the values in column 
9 assuming LTE; Column 11: Ratio of the gas mass to the virial mass; Column 12: Mass of the M01 SCUBA cores associated with each clump. 



Clump 


Size 


Peak 


0"1D 




Axis 


'Vir 




Clump 


M m 


Bound 


M smm 


No. 


/(AU X AU) 


/K 


/km s 


/kms-' 


ratio 


/pc 


/Mo 


sum /K 


/Mo 


ratio 


Ms 


1 


15000 x 13000 


10.68 


0.34 


11.04 


0.88 


0.083 


1.29 


3424 


11.89 


9.2 


2.45 


2 


14000 x 16000 


8.13 


0.36 


10.88 


0.88 


0.088 


1.27 


2960 


10.28 


8.1 


2.00 


3 


13000 x 15000 


7.06 


0.20 


10.90 


0.87 


0.082 


0.60 


1793 


6.22 


10.4 




4 


7800 x 4000 


6.59 


0.22 


10.71 


0.74 


0.038 


0.33 


546 


1.90 


5.8 




5 


17000 x 13000 


6.14 


0.42 


10.75 


0.78 


0.087 


2.72 


2645 


9.18 


3.4 


3.55 


6 


17000 x 12000 


5.78 


0.30 


10.68 


0.74 


0.084 


2.36 


1934 


6.71 


2.8 


0.60 


7 


7900 x 10000 


4.95 


0.34 


10.46 


0.82 


0.055 


1.96 


1081 


3.75 


1.9 


0.65 


8 


12000 x 9800 


5.19 


0.20 


11.01 


0.84 


0.065 


0.57 


388 


1.35 


2.4 




9 


6200 x 7700 


5.05 


0.40 


10.17 


0.88 


0.043 


0.67 


930 


3.23 


4.8 


1.30 


10 


8000 x 1 1000 


4.70 


0.21 


10.65 


0.78 


0.058 


0.75 


258 


0.89 


1.2 




1 1 


1 1000 x 9900 


4.80 


0.20 


10.42 


0.93 


0.062 


0.40 


570 


1.98 


4.9 




12 


7200 x 9700 


4.61 


0.25 


9.79 


0.83 


0.051 


0.66 


995 


3.45 


5.3 


1.90 


13 


1 1000 x 2900 


5.77 


0.27 


10.62 


0.53 


0.043 


0.61 


577 


2.00 


3.3 




14 


21000 x 14000 


4.31 


0.25 


11.02 


0.72 


0.101 


1.02 


2001 


6.95 


6.8 


1.20 


15 


9600 x 8900 


7.28 


0.31 


11.01 


0.95 


0.056 


0.47 


1257 


4.36 


9.3 


2.00 


16 


1 1000 x 7900 


3.98 


0.19 


10.91 


0.77 


0.057 


0.46 


583 


2.02 


4.4 




17 


7600 x 5700 


3.72 


0.23 


9.69 


0.85 


0.042 


0.67 


448 


1.56 


2.3 


4.30 


IS 


30000x5100 


2.82 


0.27 


10.81 


0.26 


0.083 


0.99 


575 


2.00 


2.0 


0.60 


19 


7300 x 13000 


2.86 


0.16 


10.76 


0.66 


0.059 


0.54 


577 


2.00 


3.7 




20 


7800 x 8100 


3.45 


0.42 


10.10 


0.97 


0.049 


1.26 


822 


2.85 


2.3 




21 


3400 x 7000 


2.69 


0.28 


10.12 


0.75 


0.035 


0.26 


392 


1.36 


5.3 


2.50 


22 


9500 x 17000 


2.44 


0.32 


10.91 


0.63 


0.076 


1.33 


1163 


4.04 


3.1 


2.30 


23 


3900 x 4400 


2.52 


0.24 


10.78 


0.97 


0.03 


0.34 


187 


0.65 


1.9 




24 


15000 x 4200 


2.52 


0.23 


11.11 


0.45 


0.055 


0.30 


486 


1.69 


5.7 




25 


7100 x 4600 


4.23 


0.39 


10.84 


0.82 


0.038 


0.17 


629 


2.18 


13.1 


1.35 


26 


7300 x 3200 


5.75 


0.30 


10.45 


0.72 


0.035 


0.39 


531 


1.84 


4.8 


0.90 
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